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ABSTRACT 

We study non-axisymmetric structure of accretion discs in Be/X-ray binaries, 
performing three dimensional Smoothed Particle Hydrodynamics simulations for a 
coplanar system with a short period and a moderate eccentricity. We find that ram 
pressure due to the phase-dependent mass transfer from the Be-star disc excites a 
one-armed, trailing spiral structure in the accretion disc around the neutron star. 
The spiral wave is transient; it is excited around the periastron passage, when the 
material is transferred from the Be disc, and is gradually damped afterwards. The 
disc changes its morphology from circular to eccentric with the development of the 
spiral wave, and then from eccentric to circular with the decay of the wave during one 
orbital period. It turns out that the inward propagation of the spiral wave significantly 
enhances the mass-accretion rate onto the neutron star. Thus, the detection of an X-ray 
luminosity peak corresponding to the peak in enhanced mass-accretion rate provides 
circumstantial evidence that an accretion disc is present in Be/X-ray binaries. 

Key words: accretion, accretion discs - hydrodynamics - methods: numerical - 
binaries: general - stars: emission-line, Be - X-rays: binaries 


1 INTRODUCTION 

About two-thirds of high-mass X-ray binaries have been 
identified as Be/X-ray binaries. These systems generally 
consist of a neutron star and a Be star with a cool (~ 10 4 A) 
equatorial disc, which is geometrically thin and nearly Kep- 
lerian. Be/X-xay binaries are distributed over a wide range 
of orbital periods (10 d < P or b < 300 d) and eccentricities 
(e<0.9). 

Most Be/X-ray binaries show only transient activity in 
X-ray emission and are termed Be/X-ray transients. Be/X- 
ray transients show periodic (Type I) outbursts, which 
are separated by the orbital period and have luminosity 
Lx = 10 36-3 'ergs -1 , and giant (Type II) outbursts of 
Lx > 10 37 erg s -1 with no orbital modulation. These out¬ 
bursts have features that strongly suggest the presence of 
an accretion disc around the neutron star. 

Recently, Hayasaki & Okazaki (2004, hereafter paper I) 
studied the accretion flow around the neutron star in a 
Be/X-ray binary with a short period (P or b = 24.3 d) and 
a moderate eccentricity (e = 0.34), using a three dimen¬ 
sional Smoothed Particle Hydrodynamics (SPH) code. They 
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found that a time-dependent accretion disc is formed around 
the neutron star. They also discussed the evolution of the 
azimuthally-averaged structure of the disc, although the disc 
is significantly eccentric in spite of the nealry Keplerian ro¬ 
tation. The non-axisymmetric structure seen in the disc is 
considered to result from a phase-dependent mass transfer 
from the Be disc. 

iTosa ( 4994 ) showed numerically that ram pressure due 
to an intergalactic gas excites the one-armed spiral structure 
in the non-se lfgrav itating, galactic U-constant disc. Then 
iKato fc Tosa 1 4994 ) showed analytically that the excitation 
of the one-armed mode is caused by the same resonant insta¬ 
bility mechanism as the tidally induced eccentric instability 
which e xplains the superhump phenomena in dwarf novae 
systems iHirose fc Osakl [l989j.[Lubow|l99ll). In the context 
of accretion disc theory, iLubow : ( 4994 ) investigated the ef¬ 
fect of time varying mass-transfer on the eccentricity of an 
accretion discs in a circular binary. However, it has never 
been explicitly studied how phase-dependent mass transfer 
affects the eccentricity of an accretion disc in an eccentric 
binary. 

In this letter, we study the non-axisymmetric structure 
of the accretion disc around the neutron star in Bc/X-ray 
binaries. In Section [5] the development of the one-armed 
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spiral structure is demonstrated with 3D SPH simulations. 
The strength of several modes excited in the disc, the orbital 
modulation of the disc radius and the phase dependence of 
the mass-accretion rate are also analyzed in this section. Our 
conclusions are presented in Section 0 


2 RAM-PRESSURE-DEFORMED ACCRETION 

DISCS 

In this section, we show that ram pressure due to the ma¬ 
terial transferred from the Be disc around periastron tem¬ 
porarily excites the one-armed spiral wave in the accretion 
disc around the neutron star in Be/X-ray binaries. 

Our simulations were performed using the same 3D SPH 
code as in paper I, which was based on a version originally 
developed by Benz (Benz 1990; Benz et al. 1990) and later 
by Bate, Bonnell & Price (1995). In order to investigate the 
effect of ram pressure on the accretion disc, we compare 
results from model 1 in paper I (hereafter, model A) with 
those from a simulation (hereafter, model B) in which the 
mass transfer from the Be disc is artificially stopped for one 
orbital period. Except for this difference, the two simulations 
have the same model parameters: The orbital period P or b is 
24.3 d, the eccentricity e is 0.34, and the Be disc is coplanar 
with the orbital plane. The inner radius of the simulation 
region n n is 3.Ox 10~ 3 a, where a is the semi-major axis of the 
binary. A polytropic equation of state with the exponent V = 
1.2 is adopted. The Shakura-Sunyaev viscosity parameter 
ass = 0.1 throughout the disc. There is also a bulk viscosity 
Vb = 5assCsH/3, where c 3 is the polytropic sound speed and 
H is the scale-hight of the disc. 

Fig. Ogives a sequence of snapshots of the accretion disc 
around the neutron star for 7 < t < 8 in model A, where 
the unit of time is the orbital period P or b- The left pan¬ 
els show the contour maps of the surface density, whereas 
the non-axisymmetric components of the surface density and 
the velocity field are shown in the right panels. Annotated 
in each left panel are the time in units of P or b, the relative 
disk thickness H/r at r = 0.05a and the mode strength Si, 
a measure of the amplitude of the one-armed spiral wave, 
details of which will be described later. It is noted from the 
figure that the one-armed, trailing spiral is excited at perias¬ 
tron and is gradually damped towards the next periastron. 

For comparison, we present the results for model B, 
in which the mass transfer is artificially turned off for 
7 < t < 8. Fig. 0 shows the surface density (the left panel) 
and the non-axisymmetric components of the surface den¬ 
sity and the velocity field (the right panel) at the time cor¬ 
responding to the middle panel of Fig.Q The format of the 
figure is the same as that of Fig. 0 It should be noted that 
the disc deformation due to the one-armed mode is not seen 
in model B. The disc is more circular and has a larger ra¬ 
dius in model B than in model A. This strongly suggests 
that the excitation of the one-armed spiral structure in the 
accretion disc is induced by ram pressure from the material 
transferred from the Be disc at periastron. 


2.1 Mode strength 

In this subsection, we analyse the orbital-phase dependence 
of the disc eccentricity by decomposing the surface den- 
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Figure 1 . Snapshots of the accretion disc for model A. The 
left panels show the surface density in a range of three orders of 
magnitude in the logarithmic scale, while the right panels show 
the non-axisymmetric components of the surface density (gray¬ 
scale plot) and the velocity field (arrows) in the linear scale. In 
the right panels, the region in gray (white) denotes the region with 
positive (negative) density enhancement. The periastron is in the 
^-direction and the disc rotates counterclockwise. Annotated in 
each left panel are the time in units of P or b, the relative disc 
thickness H/r at r = 0.05a and the mode strength Si. 


sity distribution into Fourier components which vary as 
exp (inuf>), wh ere m is the azimuthal harmonic number. 

Following iLubowi t l99lh . we define the mode strength 
by 

Sf,m(t) - Md ^ 2 +5rn ^ j Q r # S (P^0/(W>), (!) 

where / is either sin or cos and Md is the total disc mass 
given by J r dr f 0 rd(j>Y>(r,(f),t). Then, the total strength of 
the mode m is defined by 

S m (t) = (Scos.m + ‘S'sin,m) 1/2 - (2) 

Note that our definition of t he mode strength slightly 
differs from that of lLubow i • 199 1:. He decomposed the sur¬ 
face density distribution into double Fourier components 
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Figure 2. Same as Fig. IT! but for model B. 



Figure 3. Evolution of several non-axisymmetric modes for 0 < 
t < 8. The soild, dashed and dotted lines denote the strengths of 
the m = 1, m = 2 and m = 3 modes, respectively. 


with the azimuthal and time harmonic numbers, whereas 
we did it with only the azimuthal harmonic number. This 
is because our purpose is to analyze the time dependence 
of the non-axisymmetric structure of the accretion disc in 
detail. Note also that the definition of the mode strength is 
not unique, so Eq. gives only a rough measure of the am¬ 
plitude of each mode. If the spiral wave is tightly winding, 
Eq. 121 will underestimate the amplitude significantly. 

Fig. [3 shows the evolution of several non-axisymmetric 
modes for model A, in which the phase-dependent mass 
transfer from the Be disc is taken into account. The solid, 
dashed and dotted lines denote the strengths of m = 1, 
m = 2 and m = 3 modes, respectively. From Fig. El we 
note that the m = 1 component dominates the other com¬ 
ponents throughout the run. In the initial developing stage 
of the disc (0 < t < 5), elliptical orbits of the gas particles 
transferred from the Be disc significantly contribute to the 
strength of the m = 1 component, Si. The contribution be¬ 
comes smaller as the disc develops. Thus, the decrease in Si 
for t < 5 corresponds to the decrease in the disc eccentricity 
during this period. The decrease in eccentricity is likely to 
result from a non-zero bulk viscosity which helps to prevent 
the spontane ous growth of the eccentricity through viscous 
overstability liQgilviefcOOlh . 

In contrast, after the disc is well developed (£> 5), 
the one-armed spiral wave gives the major contribution to 
Si. The mode strength significantly modulates with orbital 
phase, corresponding to the rapid excitation and the subse¬ 
quent gradual damping of the wave. 

In order to see the time dependence of non- 
axisymmetric structure in the developed accretion disc more 
closely, we show the time dependence of the mode strength 
Si for 7 < t < 8 in the top panel of Fig. El The thick and 
thin lines denote Si in model A and in model B, respectively. 
Recall that in model B the mass transfer from the Be disc 
is artificially turned off for this period. 

In model A, the mode strength Si is dominated by the 
one-armed spiral mode. It rapidly increases to the maximum 
value of ~ 0.3 at t = 7.14 (see the middle panel of Fig.0 and 
then gradually decreases to a minimum value of ~ 0.01 at 
t = 7.8 (the bottom panel of Fig.Q. It is important to note 
that the mode strength Si in model A is much stronger and 
has much larger orbital modulation than in model B. This 
clearly shows that the one-armed mode is not excited by the 
tidal interaction with the Be star but by ram pressure of the 
material transferred from the Be disc at periastron passage. 


The middle panel of Fig. E] shows the azimuth w p be¬ 
tween the position vector of the maximum value of the non- 
axisymmetric component of the surface density and the ec¬ 
centric vector of the binary orbit (hereafter, the azimuth of 
the maximum density perturbation). As in the other panels, 
the thick and thin lines are for model A and model B, re¬ 
spectively. It should be noted that w p in both models shifts 
little over one orbital period. This is a typical characteris¬ 
tic of the one-armed modes in nearly Keplerian discs ( Kato 

EH. 

2.2 Orbital modulation in the disc radius 

In paper I, we found that the disc size modulates with the 
orbital phase; the disc shrinks at periastron and gradually 
restores its radius afterwards. Our interpretation of this or¬ 
bital modulation was that the decrease in the disc radius 
at periastron is due to a negative torque by the Be star, 
while the later expansion is due to the viscous diffusion. 
For the reason described below, this interpretation must be 
corrected. 

The bottom panel of Fig.EJshows the time dependence 
of the disc radius for 7 < t < 8, which is defined in the 
same way as in paper I (equation 5). As in the other pan¬ 
els, the thick and thin lines are for model A and model B, 
respectively. From the figure, we note that the disc radius 
in model B monotonically increases with time and shows no 
orbital modulation. This is in sharp contrast to the orbital 
modulation in the disc radius seen in model A. Since the 
only difference between model A and model B is that the 
mass transfer from the Be disc is included in model A, while 
it is not in model B, it is most likely that the reduction in 
the disc radius after periastron passage is related to the mass 
transfer from the Be disc. Apparently, the negative torque 
by the Be star has little effect on the accretion disc radius 
in these models. 

How can the mass transfer from the Be star reduce the 
accretion disc size? As shown in Fig. 3 of paper I, the circu¬ 
larization radius of the material transferred from the Be disc 
is smallest around the periastron. It is much smaller than the 
accretion disc radius shown in the bottom panel of Fig.E] In 
other words, the specific angular momentum of the material 
transferred from the Be disc is much smaller than that of the 
material at the outer radius of the accretion disc. Accretion 
of such material onto the disc outer radius will temporarily 
reduce the disc radius. This mechanism, together with the 
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Figure 4. Time dependence of the mode strength Si, the azimuth 
of the maximum density perturbation, ujp , at r = 0.05a, and the 
disc radius r\| for 7 < t < 8. In each panel, the thick line is for 
model A, while the thin line for model B. 

viscous diffusion, naturally explains the orbital modulation 
in the disc radius found in paper I. 

2.3 Phase dependence of the accretion rate 

After the accretion disc is developed (t> 5), the mass- 
accretion rate has double peaks per orbit, a relatively- 
narrow, low peak at periastron and a broad, high peak after¬ 
wards [see Fig. 15(a) of paper 1], While the first low peak at 
periastron could be artificial, being related to the presence 
of the inner simulation boundary, the origin of the second 
high peak was not clear. Below we show that the one-armed 


spiral wave is responsible for the second peak in the mass- 
accretion rate. 

Fig . 0shows the time dependence of the mass-accretion 
rate for 7 < t < 8. The thick line denotes the mass-accretion 
rate in model A, in which the mass transfer from the Be disc 
is taken into account. For comparison, the mass-accretion 
rate in model B, in which the mass transfer from the Be 
disc is artificially turned off at t = 7, is also shown by 
the thin line. The difference between the accretion rate pro¬ 
files for these two models is striking. The accretion rate in 
model B monotonically decreases over one orbital period, 
whereas that of model A shows a broad peak centred at 
t ~ 7.32 — 7.35, which corresponds to the second peak found 
in paper 1. 

Although it is obvious that the above peak is caused by 
the mass transfer from the Be disc, the mass transfer rate has 
a narrow peak at periastron. What mechanisms connecting 
the mass transfer from the Be disc and the mass accretion 
onto the neutron star lead to a phase delay of ~ 0.3? There 
are two possibilites. Either the disc is the viscously adjusting 
to the addition of low-angular momentum material. How¬ 
ever, even if the angular momentum of the material trans¬ 
ferred from the Be disc is the lowest, its viscous time scale at 
the circularization radius is over 7 binary periods [see Fig. 3 
of paper I]. Thus, the viscous accretion of the lowest-angular 
momentum material could play no key role of the enhance¬ 
ment of the accretion rate at an orbital phase of about 0.35. 
Or, the second possibility is that the one-armed spiral wave 
is excited by ram pressure of the transferred material from 
the Be disc and then the spiral wave travels inwardly. The 
excitation of the spiral wave takes about ~ (0.1 — 0.2)P or b 
(see the top panel of Fig. 0 and the difference of the peak 
position on the orbital phase between the mode strength and 
the accretion rate results from the wave propagation from 
the disc outer radius to the inner simulation boundary. 

The time-scale for wave propagation is roughly esti¬ 
mated by using the dispersion relation of th e on e-armed 
oscillation in the nearly Keplerian discs llKatoll983ft . Its fre¬ 
quency c o is written by — Q OI h{k T H) 2 /2, where fc r is the 
radial wavenumber and f2 or b is the orbital angular velocity. 
Finally, the time-scale for a global, one-armed pertubation 
to travel across the disc is 2-it/uj ~ 0.13P or b at r ~ 0.05a 
at t ~ 7.1. Therefore, combining the time-scales for wave 
excitation and propagation, we expect that the peak of the 
mass-accretion rate lags behind that of the mass-transfer 
rate by (0.2 — 0.3)P or b, which is in good agreement with the 
phase of the second peak of mass-accretion rate shown in 

Fig- El 


3 CONCLUSIONS 

In this letter, we have investigated non-axisymmetric struc¬ 
ture in the accretion disc around the neutron star in Be/X- 
ray binari es, an alysin g the res ults from a simulation per¬ 
formed bv IHavasaki fc Okazaki (1 2004 1. We have adopted the 
phase-depend ent, mass-trans f er rate from a high-resolution 
simulation by lOkazaki et al. ()200 2) for a coplanar system 
with a short period (P or b = 24.3 d) and a moderate eccen¬ 
tricity (e = 0.34), which targeted 4U0115+63, one of the 
best studied Be/X-ray binaries. 

We have found that a one-armed, trailing spiral wave 
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Figure 5. Time dependence of the mass-accretion rate for 7 < 
t < 8. The thick and thin lines are for model A and model B, 
respectively. The right axis shows the X-ray luminocity corre¬ 
sponding to the mass-accretion rate. 


is excited in the accretion disc by ram pressure from the 
strongly phase-dependent mass-transfer from the Be disc. 
The spiral wave is transient; it is excited at periastron pas¬ 
sage, when mass transfer from the Be disc takes place, and is 
gradually damped afterwards. The one-armed perturbation 
pattern in the accretion disc precesses little over one orbital 
period, which is consistent with the theory of global, one- 
armed oscillations in nearly Keplerian discs (lKat,olll983l) . 

The eccentric (m = 1) mode dominates the other non- 
axisymmetric modes throughout the run (0 < t < 8P or b). 
In the initial stage of disc development (0 < t < 5P 0 rb), 
elliptical orbits of the gas particles transferred from the Be 
disc make the accretion disc significantly eccentric. After 
the disc is well developed (t> 5), however, the one-armed 
spiral wave gives the major contribution to the eccentric 
deformation of the accretion disc. Then, the strength of the 
deformation significantly modulates with the orbital phase, 
corresponding to the rapid excitation and the subsequent 
gradual damping of the one-armed wave. 

_ The system we studied differs from that of lKato fc Tosal 

f 19941 1 in the sense that our accretion disc is nearly Keple¬ 
rian and ram pressure of the material transferred from the 
Be disc exerts on the accretion disc only briefly around peri¬ 
astron. Nevertheless, we can estimate the growth rate of the 
one-ar med spir a l wave bas ed on their formalism. Accord¬ 
ing to iKato fc Tosal f 1994h . the growth rate A is given by 
A ~ i^v r /{rA — r in) where v r is the radial velocity of the disc. 
In our case, adopting v r ~ 0.07afl or b at r ~ ra ~ 0.05a, we 
have A ~ Hk 2> fi 0 rb at r = 0.05a for 7.0 < t < 7.2. This 
explains the rapid growth of the spiral wave for 7.0 <t< 7.2 
in our simulation. 

After the accretion disc has developed (t>5P or b), the 
mass-accretion rate has double peaks per orbit, a relatively- 
narrow, low peak at periastron and a broad, high peak af¬ 
terwards. We have found that the one-armed spiral wave 
excited by ram pressure of the material transferred from 


the Be disc enhances the inward mass flux in the accretion 
disc, and is therefore responsible for this high peak in the 
accretion rate. The X-ray luminosity peak corresponding to 
the high peak in accretion rate has never been identified 
in Be/X-ray binaries. However, if this peak is detected, it 
could be understood as strong circumstantial evidence for 
the presence of an accretion disc in Be/X-ray binaries. 

In this letter, we have reported on the effects of ram 
pressure on the accretion disc, including the excitation of 
the one-armed spiral wave, which enhances the accretion 
rate. These results are based on the simulations, of which 
the run time was much shorter than the viscous time-scale. 
In longer simulation runs, new features/interactions are ex¬ 
pected to appear. While the deformation of the accretion 
disc by ram pressure is expected to become small as the disc 
mass increases, the excitation of the eccentric mode directly 
driven by the one-armed bar potential, which is not seen in 
our simulations, will become important. The tidal/resonant 
torque by the Be star may begin to work on the outer part 
of the accretion disc, as the disc size increases. In a subse¬ 
quent paper, we will investigate the long-term evolution of 
the accretion disc in Be/X-ray binaries in detail. 
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